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We present an updated measurement of the anomalous like-sign dimuon charge asymmetry 
for semi-leptonic 6-hadron decays in 9.0 fb~^ of pp collisions recorded with the DO detector at 
a center-of-mass energy of = 1.96 TeV at the Fermilab Tevatron collider. We obtain = 
(-0.787 ± 0.172 (stat) ± 0.093 (syst))%. This result differs by 3.9 standard deviations from the 
prediction of the standard model and provides evidence for anomalously large CP violation in semi- 
leptonic neutral B decay. The dependence of the asymmetry on the muon impact parameter is 
consistent with the hypothesis that it originates from semi-leptonic 6-hadron decays. 

PACS numbers: 13.25.Hw; 14.40.Nd; ll.30.Er 



I. INTRODUCTION 

We measure the like-sign dimuon charge asymmetry of 
semi-leptonic decays of b hadrons, 

in 9.0 fb~^ oipp collisions recorded with the DO detector 
at a center-of-mass energy ^/s = 1.96 TeV at the Fermi- 
lab Tevatron collider. Here N^^ and arc the num- 
ber of events containing two positively charged or two 



*with visitors from "Augustana College, Sioux Falls, SD, USA, 
''The University of Liverpool, Liverpool, UK, '^SLAC, Mcnlo Park, 
CA, USA, '^University College London, London, UK, '^Ccntro 
de Invcstigacion en Computacion - IPN, Mexico City, Mexico, 
•'^ECFM, Universidad Autonoma de Sinaloa, Culiacan, Mexico, and 
^Universitiit Bern, Bern, Switzerland. ^Deceased. 



negatively charged muons, respectively, both of which are 
produced in prompt semi-leptonic 6-hadron decays. At 
the Fermilab Tevatron pp collider, h quarks are produced 
mainly in hh pairs. Hence, to observe an event with two 
like-sign muons from semi-leptonic 5-hadron decay, one 
of the hadrons must be a i?° or B° meson that oscillates 
and decays to a muon of charge opposite of that expected 
from the original b quark [1]. The oscillation O ZJj^ 
{q = d ov s) is described by higher order loop diagrams 
that are sensitive to hypothetical particles that may not 
be directly accessible at the Tevatron. 

The asymmetry ^1^^ has contributions from the semi- 
leptonic charge asymmetries a^j and of and B^ 
mesons [2], respectively: 

Al, = Caai + Csal,, (2) 

Ar 

witha^i = ^^tan^q, (3) 
where 0^ is a CP-violating phase, and AMg and AF^ are 
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the mass and width differences between the eigenstates of 
the propagation matrices of the neutral mesons. The 
coefficients Cd and Cs depend on the mean mixing proba- 
bility, Xo, and the production rates of B° and B° mesons. 
We use the values of these quantities measured at LEP as 
averaged by the Heavy Flavor Averaging Group (HFAG) 
[3] and obtain 



Cd 



0.594 ±0.022, 
0.406 ±0.022. 



(4) 



The value of xo measured by the CDF Collaboration re- 
cently [4] is consistent with the LEP value, which sup- 
ports this choice of parameters. Using the standard 
model (SM) prediction for a^i and afj [5], we find 



Al,{SM) = (-0.028troo6) 



(5) 



which is negligible compared to present experimental sen- 
sitivity. Additional contributions to CP violation via 
loop diagrams appear in some extensions of the SM and 
can result in an asymmetry A^^ within experimental reach 
[6-10]. 

This Article is an update to Ref. [11] that reported ev- 
idence for an anomalous like-sign dimuon charge asym- 
metry with 6.1 fb^^ of data, at the 3.2 standard devia- 
tion level. All notations used here are given in Ref. [11]. 
This new measurement is based on a larger dataset and 
further improvements in the measurement technique. In 
addition, the asymmetry's dependence on the muon im- 
pact parameter (IP) [12] is studied. The DO detector 
is described in Ref. [13]. We include a brief overview 
of the analysis in Sec. II. Improvements made to muon 
selections are presented in Sec. Ill; the measurement of 
all quantities required to determine the asymmetry A^^ 
is described in Sees. IV X. and the result is given in 
Sec. XI. Sections XII XIII present consistency checks of 
the measurement; Sec. XIV describes the study of the 
asymmetry's IP dependence. Conclusions are given in 
Sec. XV. 



II. METHOD 

The elements of our analysis are described in detail in 
Ref. [11]. Here, we summarize briefly the method, em- 
phasizing the improvements to our previous procedure. 
We use two sets of data; (i) inclusive muon data collected 
with inclusive muon triggers that provide positively 
charged muons and n~ negatively charged muons, and 
(ii) like-sign dimuon data, collected with dimuon trig- 
gers that provide N~^~^ events with two positively charged 
muons and A'' events with two negatively charged 
muons. If an event contains more than one muon, each 
muon is included in the inclusive muon sample. Such 
events constitute about 0.5% of the total inclusive muon 
sample. If an event contains more than two muons, the 
two muons with the highest transverse momentum (pt) 
are selected for inclusion in the dimuon sample. Such 



events comprise about 0.7% of the total like-sign dimuon 
sample. 

From these data wo obtain the inclusive muon charge 
asymmetry a and the like-sign dimuon charge asymmetry 
A, defined as 



a = 



A = 



n+ + n~ 



N- 



N++ + N- 



(6) 



In addition to a possible signal asymmetry ^^j, these 
asymmetries have contributions from muons produced in 
kaon and pion decay, or from hadrons that punch through 
the calorimeter and iron toroids to penetrate the outer 
muon detector. The charge asymmetry related to muon 
detection and identification also contributes to a and A. 
These contributions are measured with data, with only 
minimal input from simulation. The largest contribution 
by far is from kaon decays. Positively charged kaons have 
smaller cross sections in the detector material than nega- 
tively charged kaons [14] , giving them more time to decay. 
This difference produces a positive charge asymmetry. 

We consider muon candidates with px in the range 1.5 
to 25 GeV. This range is divided into six bins as shown 
in Table I. The inclusive muon charge asymmetry a can 
be expressed [11] as 

6 

a = E + S,) + fkai, + flai + (7) 

i=l 

where the fraction of reconstructed muons, in a given 
Pt interval i in the inclusive muon sample is given in Ta- 
ble I. The fractions of these muons produced by kaons, 
pions, and protons in a given px interval i are f]^, fl, 
and fp, and their charge asymmetries are a}^, a]^, and 
Up, respectively. We refer to these muons as "long" or 
"i" muons since they are produced by particles travel- 
ing long distances before decaying within the detector 
material. The track of a L muon in the central tracker is 
dominantly produced by the parent hadron. The charge 
asymmetry of L muons results from the difference in the 
interactions of positively and negatively charged parti- 
cles with the detector material, and is not related to 
CP violation. The background fraction is defined as 
/^kg = fk + n + Pp- The quantity /| = 1 - /^^g is the 
fraction of muons from weak decays of h and c quarks 
and r Icptons, and from decays of short-lived mesons 
{(j),ui,r], p°). We refer to these muons as "short" or "S"' 
muons, since they arise from the decay of particles at 
small distances from the pp interaction point. These par- 
ticles are not affected by interactions in the detector ma- 
terial, and once muon detection and identification imbal- 
ances are removed, the muon charge asymmetry as must 
therefore be produced only through CP violation in the 
underlying physical processes. The quantity Si in Eq. (7) 
is the charge asymmetry related to muon detection and 
identification. The background charge asymmetries a^. 
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a^, and arc measured in tlie inclusive muon data, and 
include any detector asymmetry. The therefore ac- 
counts only for S muons and is multiplied by the factor 
fh- 

The like-sign dimuon charge asymmetry A can be ex- 
pressed [11] as 

6 

A = FssAs + FsLas + ^ f;{{2 - F{,^^)Si 

i=l 

+F^j,ai, + Flai + F;ai,}. (8) 

The quantity As is the charge asymmetry of the events 

with two like-sign S muons. The quantity ^55 is the 
fraction of like-sign dimuon events with two S muons, 
FsL is the fraction of like-sign dimuon events with one 
S and one L muon. Wc also define the quantity Fj^i^ 
as the fraction of like-sign dimuon events with two L 
muons. The quantity i^^ is the fraction of muons in the 
Pt interval i in the like-sign dimuon data. The quantities 
(x = K, Tr,p) are defined as = 2Af*/iV;, where iV^ 
is the number of muons produced by kaons, pions, and 
protons, respectively, in apx interval i, with NJ^ being the 
number of muons in this interval, with the factor of two 
taking into account the normalization of these quantities 
per like-sign dimuon event. The quantity ^^kg ^ sum 
over muons produced by hadrons: 

i^bkg ^ + + F;. (9) 
We also define Fbkg as 

6 

Fbkg = E(^M^bkg) (10) 

i=l 

= Fsl + 2Fll 

= 1 + Fll-Fss. (11) 

The estimated contribution from the neglected quadratic 

terms in Eq. (8) is approximately 2 x 10^"'', which cor- 
responds to about 5% of the statistical uncertainty on 
A. 

The asymmetries as and in Eqs. (7) and (8) are 
the only asymmetries due to CP violation in the pro- 
cesses producing S muons, and are proportional to the 
asymmetry A^^: 

as = CfcAgi, 

As = CbAl,. (12) 

The dilution coefficients Ch and C(, are discussed in 
Ref. [11] and in Sec. X below. 

Equations (7) - (12) are used to measure the asym- 
metry Agj. The major contributions to the uncertainties 
on are from the statistical uncertainty on A and the 
total uncertainty on Fj^, f]^ and J,. To reduce the latter 
contributions, we measure the asymmetry ^4^^ using the 
asymmetry A' , which is defined as 

A' = A - aa. (13) 



TABLE I: Fractions of muon candidates in the inclusive muon 
sample (/^) and in the like-sign dimuon sample (F^, with two 
entries per event). 



Bin Muon pr range (GeV) fi 

1 1.5 - 2.5 0.0077 0.0774 

2 2.5 - 4.2 0.2300 0.3227 

3 4.2 - 5.6 0.4390 0.3074 

4 5.6 - 7.0 0.1702 0.1419 

5 7.0 - 10.0 0.1047 0.1057 

6 10.0 - 25.0 0.0484 0.0449 



Since the same physical processes contribute to both F]^ 
and their uncertainties are strongly correlated, and 
therefore partially cancel in Eq. (13) for an appropriate 
choice of the coefficient a. The contribution from the 
asymmetry A^j, however, does not cancel in Eq. (13) be- 
cause Cft ^ Cb [11]. Pull details of the measurements of 
different quantities entering in Eqs. (7) (12) arc given 
in Ref. [11]. The main improvements in the present anal- 
ysis are related to muon selection and the measurement 
of and /If. These modifications are described in Sec- 
tions 111, IV and V. 



III. MUON SELECTION 

The muon selection is similar to that described in 
Ref. [11]. The inclusive muon and like-sign dimuon sam- 
ples are obtained from data collected with single and 
dimuon triggers, respectively. Charged particles with 
transverse momentum in the range 1.5 < px < 25 GeV 
and with pseudorapidity |r?| < 2.2 [15] are considered as 
muon candidates. The upper limit on pT is applied to 
suppress the contribution of muons from W and Z bo- 
son decays. To ensure that the muon candidate passes 
through the detector, including all three layers of the 
muon system, we require either pt > 4.2 GeV or a lon- 
gitudinal momentum (;omponcnt \pz\ > 5.4 GeV. Muon 
candidates are selected by matching central tracks with 
a segment reconstructed in the muon system and by ap- 
plying tight quality requirements aimed at reducing false 
matching and background from cosmic rays and beam 
halo. The transverse impact parameter of the muon track 
relative to the reconstructed pp interaction vertex must 
be smaller than 0.3 cm, with the longitudinal distance 
from the point of closest approach to this vertex smaller 
than 0.5 cm. Strict quality requirements are also applied 
to the tracks and to the reconstructed pp interaction ver- 
tex. The; inclusive muon sample contains all muons pass- 
ing the selection requirements. If an event contains more 
than one muon, each muon is included in the inclusive 
muon sample. The like-sign dimuon sample contains all 
events with at least two muon candidates with the same 
charge. These two muons are required to have an invari- 
ant mass greater than 2.8 GeV to minimize the number 
of events in which both muons originate from the same b 
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quark (e.g., 6 — >• /z, & — >■ c — >■ /z). Compared to Ref. [11], 
the following modifications to the muon selection are ap- 
plied: 

• To reduce background from a mismatch of tracks 
in the central detector with segments in the outer 
muon system, we require that the sign of the cur- 
vature of the track measured in the central tracker 
be the same as in the muon system. This selec- 
tion was not applied in Ref. [11], and removes only 
about 1% of the dimuon events. 

• To ensure that the muon candidate can penetrate 
all three layers of the muon detector, we require 
either a transverse momentum > 4.2 GeV, 
or a longitudinal momentum component \pz\ > 
5.4 GeV, instead of pt > 4.2 GeV or > 6.4 
GeV in Ref. [11]. With this change, the number 
of like-sign dimuon events increases by 25%, with- 
out impacting the condition that the muon must 
penetrate the calorimeter and toroids, as can be 
deduced from Fig. 1. 

• To reduce background from kaon and pion decays 
in flight, we require that the calculated from the 
difference between the track parameters measured 
in the central tracker and in the muon system be 

< 12 (for 4 d.o.f.) instead of 40 used in Ref. [11]. 
With this tighter selection, the number of like-sign 
dimuon events is decreased by 12%. 

Compared to the selections applied in Ref. [11], the 
total number of like-sign dimuon events after applying 
all these modifications is increased by 13% in addition 
to the increase due to the larger integrated luminosity of 
this analysis. 

The muon charge is determined by the central tracker. 
The probability of charge mis-measurement is obtained 
by comparing the charge measured by the central tracker 
and by the muon system and is found to be less than 
0.1%. 

The polarities of the toroidal and solenoidal magnetic 
fields are reversed on average every two weeks so that the 
four solenoid-toroid polarity combinations are exposed to 
approximately the same integrated luminosity. This al- 
lows for a cancellation of first-order effects related to the 
instrumental asymmetry [16]. To ensure such cancella- 
tion, the events are weighted according to the number 
of events for each data sample corresponding to a dif- 
ferent configuration of the magnets' polarities. These 
weights are given in Table II. During the data taking 
of the last part of the sample, corresponding to approx- 
imately 2.9 fb~^ of pp collisions, the magnet polarities 
were specially chosen to equalize the number of dimuon 
events with different polarities in the entire sample. The 
weights in Table II are therefore closer to unity compared 
to those used in Ref. [11]. 




P, [GeV] 



FIG. 1: (color online). Smallest muon momentum required 
to penetrate the calorimeter and toroids at different pseudo- 
rapidities, \ri\ (solid line), and the momentum selection used 
in this analysis (dashed line). 



TABLE II: Weights assigned to the events recorded with dif- 
ferent solenoid and toroid polarities in the inclusive muon and 
like-sign dimuon samples. 



Solenoid 


Toroid 


Weight 


Weight 


polarity 


polarity 


inclusive muon 


like-sign dimuon 


-1 


-1 


0.994 


0.964 


-1 


-1-1 


1.000 


1.000 


4-1 


-1 


0.985 


0.958 


-1-1 


+ 1 


0.989 


0.978 



IV. MEASUREMENT OF /jc, U, fp 

The fraction /J^- in the inclusive muon sample is mea- 
sured using K*'-' K'^TT~ decays, with the kaon identi- 
fied as a muon (see Ref. [11] for details). The transverse 
momentum of the meson is required to be in the pt 
interval i. Since the momentum of a particle is measured 
by the central tracking detector, a muon produced by a 
kaon is assigned the momentum of this kaon (a small cor- 
rection for kaons decaying within the tracker volume is 
introduced later). The fraction /]^.o of these decays is 
converted to the fraction using the relation 



K 



where ni{Kg) and ni{K*^ 



(14) 



reconstructed Kg 



KgTT^) are the number of 



TT and K*+ — >• K*^tt^ decays, 
The transverse momentum of the Kg me- 



respectively. 
son is required to be in the px interval i 



We require in 
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FIG. 2: (color online). The tt^tt" invariant mass distribu- 
tion for Kg candidates in the inclusive muon sample with at 
least one pion identified as a muon with 4.2 < Pt{K'^) < 5.6 
GeV. The solid line represents the result of the fit to the 
content, and the dashed line represents the fitted background 
contribution. 



addition that one of the pions from the Kg — ;> tt+tt" 
decay be identified as a muon. In the previous analy- 
sis [11] the production of K*^ mesons was studied in a 
sample of events with an additional reconstructed muon, 
but we did not require that this muon be associated with 
a pion from Kg tt+tt^ decay. The fraction of events 
containing b and/or c quarks was therefore enhanced in 
the sample, which could result in a bias of the measured 
fraction Jk. This bias does not exceed the systematic un- 
certainty of Jk and its impact on the A^^ value is less than 
0.03%. The application of the new requirement ensures 
that the flavor composition in the selected K*'^ — Kgir'^ 
and K*'^ — > K'^tt~ samples is the same and this bias is 
eliminated. 

The selection criteria and fitting procedures used to 



select and determine the number of Kg, K*'^ and K*^ 
events are given in Rcf. [11]. As an example. Fig. 2 
displays the tt+tt^ invariant mass distribution and the 
fitted Kg — ^ TT+TT^ candidates in the inclusive muon 
sample, with at least one pion identified as a muon, for 
4.2 < pt{K^s) < 5-6 GeV. Figure 3 shows the K'^n+ 
mass distribution and fit to K*~^ — )• Kgir'^ candidates 
for all /eg candidates with 4.2 < priK^) < 5.6 GeV and 
480 < M(7r+7r-) < 515 MeV. Figure 4 shows the K+ir' 
mass distribution and the fit result for K*'^ — > A'+tt^ 
candidates for all kaons with 4.2 < pt{K~^) < 5.6 GeV. 
The K~^TT~ mass distribution contains contributions from 
light meson resonances decaying to tt+tt^. The most 
important contribution comes from the p'^ — > tt+tt" de- 
cay with TT — > /i. It produces a broad peak in the mass 
region close to the K*'^ mass. The distortions in the 
background distribution due to other light resonances, 
which arc not identified explicitly, can also be seen in 
Fig. 4. Our background model therefore includes the 
contribution of p'^ — >■ 7r"'"7r~ and two additional Gaus- 
sian terms to take into account the distortions around 
1.1 GeV. More details of the background description are 
given in Ref. [11]. 
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FIG. 3: (color online), (a) The K'^tt'^ invariant mass dis- 
tribution for K*'^ candidates in the inclusive muon sample. 
The Kg candidate is required to have 480 < M(7r"'"7r~) < 515 
MeV and 4.2 < Pt{Ks) < 5.6 GeV. The solid line represents 
the result of the fit to the K*^ content, and the dashed line 
shows the background contribution, (b) Difference between 
data and the result of the fit. 



The measurement of the fractions /^r and fp is also 
performed using the method of Ref. [11]. The values of 
fx and /tt are divided by the factors Ck and C^r , respec- 
tively, which take into account the fraction of kaons and 
pions reconstructed by the tracking system before they 
decay. These factors are discussed in Ref. [11], and are 
determined through simulation. Contrary to Ref. [11], 
this analysis determines these factors separately for kaons 
and pions. We find the values: 



Ck ^ 0.920 ±0.006, 
= 0.932 ±0.006. 



(15) 



The uncertainties include contributions from the num- 
ber of simulated events and from the uncertainties in the 
momentum spectrum of the generated particles. 

The values of /k , fir and fp in different muon px bins 
are shown in Fig. 5 and in Table III. The changes in 
the muon candidates selection adopted here is the main 
source of differences relative to the corresponding values 
in Ref. [11]. The fractions /,r and fp are poorly measured 
in bins 1 and 2, and bins 5 and 6 due to the small num- 
ber of events, and their contents are therefore combined 
through their weighted average. 
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FIG. 4: (color online), (a) The K'^n invariant mass dis- 
tribution for K*" candidates in the inclusive muon sample 
for all kaons with 4.2 < pt{K+) < 5.6 GeV. The solid line 
corresponds to the result of the fit to the content, and 
the dashed line shows the contribution from combinatorial 
background. The shaded histogram is the contribution from 
p^' — >■ tt'^it" events, (b) Difference between data and the re- 
sult of the fit. 



TABLE III: Fractions /k, fir, and fp for different pt bins. 
The bottom row shows the weighted average of these quan- 
tities obtained with weights given by the fraction of muons 
in a given pr interval, in the inclusive muon sample, see 
Table I. Only statistical uncertainties are given. 



Bin 


U X 10^ 


U X 10' 


fp X 10' 


1 
2 


9.35 ± 4.77 
14.91 ± 1.00 


36.20 ±4.12 


0.55 ±0.24 


3 


16.65 ± 0.41 


31.42 ±2.57 


0.11 ±0.29 


4 


17.60 ± 0.49 


27.41 ±3.46 


0.63 ±0.58 


5 
6 


14.43 ± 0.45 
12.75 ± 0.97 


19.25 ±3.19 


0.64 ±0.71 


All 


15.96 ± 0.24 


30.01 ± 1.60 


0.38 ±0.17 



V. MEASUREMENT OF Fk, F„, Fp 



The quantity Fk is expressed as 



Fk = RkJk, 



(16) 



where Rk is the ratio of the fractions of muons produced 
by kaons in like-sign dimuon and in inclusive muon data. 
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FIG. 5: (color online). The fraction of (a,) K ^ jj, tracks, (b) 
n fi tracks and (c) p — >■ /i tracks in the inclusive muon sam- 
ple as a function of the kaon, pion and proton pr, respectively. 
The horizontal dashed lines show the mean values. 



For the pt interval i, Rk is defined as 

^ N^jK^fi) mill) 

where iV, (K — >■ /i) and n,; {K — >■ /i) are the number of re- 
constructed K mesons identified as muons in the like-sign 
dimuon and in the inclusive muon samples, respectively. 
The transverse momentum of the K meson is required to 
be in the pr interval i. The quantities Ni{p) and ni{p) 
are the number of muons in the pT interval i. A mul- 
tiplicative factor of two is included in Eq. (17) because 
there are two muons in a like-sign dimuon event, and Fk 
is normalized to the number of like-sign dimuon events. 

In the previous analysis [11], the quantity Fk was ob- 
tained from a measurement of the K*'^ production rate. 
Presenting it in the form of Eq. (16) also allows the de- 
termination of Fk through an independent measurement 
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of the fraction of Kg mesons in dimuon and in inclusive 
muon data where one of the pions from Kg — >■ '!t~^tt~ 
decay is identified as a muon. This measurement is dis- 
cussed below. In addition, Eq. (16) offers an explicit sep- 
aration of systematic uncertainties associated with Fk- 
The systematic uncertainty on the fraction /k affects 
the two determinations of A^^ based on Eqs. (7) and (8) 
in a fully correlated way; therefore, its impact on the 
measurement obtained using Eq. (13) is significantly re- 
duced. The systematic uncertainty on the ratio Rk does 
not cancel in Eq. (13). It is estimated directly from a 
comparison of the values of Rk obtained in two indepen- 
dent channels. 

One way to measure Rk is from the fraction of K*^ — > 
K~^t:~ events in the inclusive muon and like-sign dimuon 
data, 



(18) 



where Ni{K*^ fi) and ni{K*^ — /i) are the number 
of reconstructed K*'^ — > X+vr^ decays, with the kaon 
identified as a muon in the like-sign dimuon and in the 
inclusive muon samples, respectively. The transverse mo- 
mentum of the K meson is required to be in the pt in- 
terval i. The measurement using Eq. (18) is based on the 
assumption 



n,{K*°^fi) n,{K^fi)' 



(19) 



which was validated through simulations in Ref. [11]. The 
corresponding systematic uncertainty is discussed below. 

In Ref. [11], the fractions Fff.o and /x*o were obtained 
independently from a fit of the K'^'!t~ invariant mass 
distribution in the like-sign dimuon and inclusive muon 
sample, respectively. Figure 6 shows the same mass stud- 
ies as in Fig. 4, but for the like-sign dimuon sample. The 
fit in both cases is complicated by the contribution from 
light meson resonances that decay to tt+tt", producing 
a reflection in the K~^Tr~ invariant mass distribution. In 
addition, the detector resolution is not known a priori 
and has to be included in the fit. All these complica- 
tions are reduced significantly or eliminated in the "null- 
fit" method introduced in Ref. [11], which is used in this 
analysis to measure the ratio Rk{K*^). 

In this method, for eachpy interval i, we define a set of 
distributions Pi{MK-K]£,) that depend on a parameter ^: 



P^{MK.■,O^N,{MK.) 



(20) 



where Ni{MKT^) and ni{MKT^) are the number of en- 
tries in the pr bin i of the K'^'n~ invariant mass dis- 
tributions in the like-sign dimuon and inclusive muon 
samples, respectively. For each value of ^ the number 
of K*^ K+Tr~ decays, J\f{K*"), and its uncertainty, 
AJ\f{K*^), are measured from the Pi{MKn',i) distribu- 
tion. The value of ^ for which J\f{K*°) = defines 
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FIG. 6: (color online), (a) The K^n invariant mass dis- 
tribution of candidates in the like-sign dimuon sample 
for all kaons with 4.2 < pt{K+) < 5.6 GeV. The solid line 
corresponds to the result of the fit to the K*'^ content, and 
the dashed line shows the contribution from combinatorial 
background. The shaded histogram is the contribution from 
p^' — > tt^tt" events, (b) Difference between data and the re- 
sult of the fit. 



RK,i{K*°). The uncertainty a{RK,i) is determined from 
the condition that N{K*^) = ±I^{K*°) corresponding 
toe = i?if,i(i^*°)±a(i?K,.). 

The advantage of this method is that the influence 
of the detector resolution becomes minimal for M{K*^) 
close to zero, and the contribution from the peaking back- 
ground is reduced in Pi {Mktt\ f o fl^*^ same extent as the 
contribution of K*^ mesons, and becomes negligible when 
J\f{K*^) is close to zero. As an example, Fig. 7 shows the 
mass distribution Pi{MK-K\0 for ^ = 0.88, for all kaons 
with 4.2 < pt{K+) < 5.6 GeV. This distribution is ob- 
tained from the distributions shown in Figs. 4 and 6, 
using Eq. (20). The contributions of both K*° K+tt' 
and TT+TT", as well as any other resonance in the 

background, disappear. As a result, the fltting procedure 
becomes more robust, the fitting range can be extended, 
and the resulting value of Rk{K*^) becomes stable under 
a variation of the fitting parameters over a wider range. 

The value of Rk is also obtained from the production 
rate of if^ mesons in the inclusive muon and dimuon 
samples. We compute RK,i for a given pT interval i, as 



''^'^ n,iK° ^ /i) N,i^i) 



(21) 



where Ni{Kg fi) and ni{Kg ji) are the num- 
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FIG. 7: (color online), (a) The K^n invariant mass distri- 
bution P2{Mk^;0 obtained using Eq. (20) for ^ = 0.88 for 
all kaons with 4.2 < pt{K'^) < 5.6 GeV. The dashed hue 
shows the contribution from the combinatorial background, 
(b) Difference between data and the result of the fit. 



TABLE IV: Values of k in different pr bins. The bottom row 
shows their average. Only statistical uncertainties are given. 

bin K 



Mean 



1.005 
1.025 
1.038 
1.036 
1.051 
1.080 



0.024 
0.016 
0.016 
0.016 
0.016 
0.013 



1.046 ± 0.007 
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FIG. 8: (color online). The correction coefficient 
function of the kaon transverse momentum. The horizontal 
dashed line shows the mean value. 



ber of reconstructed A'^ — > tt+tt^ decays with one pion 
identified as a muon in the dimuon and the inclusive 
muon data, respectively. The correction factor Ki is dis- 
cussed later in this section. The measurement of RK,i 
using Eq. (21) assumes isospin invariance and consequent 
equality of the ratio of production rates in the dimuon 
and in the inclusive muon samples of and Kg mesons, 
i.e.. 

Since the charged kaon px in Eq. (22) is required to be 
within the pt interval i, the transverse momentum of the 
Kg meson in Eq. (21) is also required to be within the px 
interval i. We expect approximately the same number of 
positive and negative pions from Kg — ^ tt+tt" decays to 
be identified as a muon. Therefore, we use both like-sign 
and opposite-sign dimuon events to measure Ni {Kg — /i) 
and we do not use the multiplicative factor of two in 
Eq. (21). The requirement of having one pion identified 
as a muon makes the flavor composition in the samples 
of charged K ^ events and Kg — >■ fi events similar. 

The charges of the kaon and the additional muon 
in a dimuon event can be correlated, i.e., in general 
N{K+^i+) ^ N{K-^i+). However, the number of 
Ni{Kg — >• fi) events is not correlated with the charge of 
the additional muon, i.e., N{K^, — /x+,/Lt+) = N{K^ — 



/i^, /i+). Since the ratio RK,i is determined for the sam- 
ple of like-sign dimuon events, we apply in Eq. (21) the 
correction factor Ki, defined as 

^ 2iNiK+fi+) + c.c.) 

'~ {N{K+fi+)+N{K-fi+) + c.cy ^ ' 

to take into account the correlation between the charges 
of the kaon and muon. The abbreviation "c.c." in 
Eq. (23) denotes "charge conjugate states". The coef- 
ficients Ki are measured in data using the events with 
a reconstructed K*^ K^n" decay and an additional 
muon. To reproduce the selection for the dimuon sample 
[11], the invariant mass of the K^ system, with the kaon 
assigned the mass of a muon, is required to be greater 
than 2.8 GeV. The fitting procedure and selection crite- 
ria to measure the number of K*^ events are described 
in Ref. [11]. The values of Ki for different pr intervals 
are given in Fig. 8 and in Table IV. 

The average muon detection efficiency is different for 
the inclusive muon and like-sign dimuon samples be- 
cause of different pt thresholds used in their triggers. 
The difference in muon detection efficiency is large for 
muons with small pT^ but it is insignificant for muons 
above the inclusive-muon trigger threshold. The ratio 
Ni{K^ — > fj,)/ni{K^ — > fi) in Eq. (21) is measured as 
a function of the transverse momenta of Kg mesons, 
pxiKg), while the ratio ni{p)/Ni{^) is measured in bins 
of muon Pt- Each px^Kg) bin contains tt ^ /i with dif- 
ferent ptIt^ — > m) values. The muon detection efficiency 
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FIG. 9: (color online). The ratio of tt — >■ detection efficien- 
cies for the inclusive muon and dimuon data as a function of 
the muon transverse momentum. The horizontal dashed line 
shows the mean value for pt{K) > 4.2 GeV. 



therefore does not cancel in Eq. (21), and can affect the 
measurement of Rk{K^). Figure 9 shows the ratio of 
TT — > /i detection efficiencies in the inclusive muon and 
dimuon data. To compute this ratio, we select the 
mesons in a given pT^K'g) interval. The pt{t^) distribu- 
tion of pions produced in the K'g — tt+tt" decay with 
a given pxiKg) is the same in the dimuon and inclusive 
muon data. Therefore, any difference in this priT^ fJ-) 
distribution between dimuon and inclusive muon data is 
due to the tt — >• detection. We compute the ratio of 
these ptIt^ — >■ fJ-) distributions, and normalize it such 
that it equals unity for pri"^ fj.) > 5.6 GeV. The value 
of this pt threshold corresponds to the px threshold for 
single muon triggers. Figure 9 presents the average of 
the ratios for different prip-) intervals. The ratio is sup- 
pressed for pt{t^ ^ /^) < 4.2 GeV, and is consistent with 
a constant for pt{t^ — >■ /i) > 4.2 GeV. To remove the 
bias due to the trigger threshold, we measure Rk{K^) 
for events with pri"^ fi) > 4.2 GeV. As a result, the 
ratio Rk is not defined for the first two pr bins in the 
Kg channel. 



The values of Rk{K*^) obtained through the null- 
fit method, for different muon pT bins, are shown in 
Fig. 10(a) and in Table V. The values of Rk{KI) are 
contained in Fig. 10(b) and in Table V. The difference 
between the values of Rk measured with K*^ mesons and 
with K'g mesons is shown in Fig. 11. The mean value of 
this difference is 



Ai?A' = 0.01 ±0.05, 



(24) 



and the x^/d.o.i. is 1.7/4. We use two independent 
methods, each relying on different assumptions, to mea- 
sure the ratio Rk and obtain results that are consistent 
with each other. The methods are subject to different 
systematic uncertainties, and therefore provide an im- 
portant cross-check. As an independent cross-check, the 
value of Rk obtained in simulation is consistent with 
that measured in data, see Sec. XIII for details. We take 
the average of the two channels weighted by their uncer- 
tainties as our final values of Rk for pt{K) > 4.2 GeV 



TABLE V: Values of Rk obtained using K*° and K% meson 
production in difi'erent pr bins. The bottom row shows their 
average. Only statistical uncertainties are given. The ratio 
Rk in the K'g channel is not measured in the first two bins, 
see Sec. V. 



bin 


Rk from K'^ 


Rk from K% 


average Rk 


1 


0.983 ± 0.154 




0.983 ± 0.154 


2 


0.931 ± 0.058 




0.931 ± 0.058 


3 


0.880 ± 0.052 


0.844 ± 0.059 


0.864 ± 0.039 


4 


0.856 ± 0.082 


0.800 ± 0.040 


0.811 ± 0.036 


5 


0.702 ± 0.112 


0.828 ± 0.042 


0.813 ± 0.039 


6 


1.160 ± 0.165 


1.138 ± 0.117 


1.146 ± 0.095 


Mean 


0.892 ± 0.032 


0.834 ± 0.025 


0.856 ± 0.020 


TABLE VI: Values of Fk, F^, and Fp for different pr bins. 


The last line shows the wei 


ghted average of these quantities 


obtained with weights given by the fraction 


of muons in a 


given Pt 


interval FJ^ in the dimuon sample, see Table I. Only 


statistical uncertainties are 


given. 




Bin 


Fk X 10^ 


F^ X 10^ 


Fp X 10^ 


1 
2 


9.19 ± 4.90 
13.88 ± 1.26 


30.54 ± 3.89 


0.47 ± 0.21 


3 


14.38 ± 0.74 


24.43 ± 2.28 


0.09 ±0.22 


4 


14.26 ± 0.74 


19.99 ±2.67 


0.46 ± 0.42 


5 
6 


11.73 ± 0.67 
14.48 ± 1.64 


14.90 ±2.55 


0.49 ±0.55 


AU 


13.78 ± 0.38 


24.81 ± 1.34 


0.35 ±0.14 



and use the values measured in the K*'^ channel for 
Pt{K) < 4.2 GeV. These values are given in Table V 
and in Fig. 10(c). As we do not observe any difference 
between the two measurements, we take half of the un- 
certainty of ARk as the systematic uncertainty of Rk- 
This corresponds to a relative uncertainty of 3.0% on the 
value of Rk- In our previous measurement [11], this un- 
certainty was 3.6%, and was based on simulation of the 
events. 

Using the extracted values of Rk, we derive the val- 
ues of Fk, and Fp. The computation of Fk is done 
using Eq. (16), and we follow the procedure described in 
Ref. [11] to determine and Fp. The results are shown 
in Fig. 12 and in Table VI. The fractions F,r and Fp are 
poorly determined for the lowest and highest pr because 
of the small number of events. The content of bins 1 and 
2, and bins 5 and 6 are therefore combined. 



VI. SYSTEMATIC UNCERTAINTIES FOR 
BACKGROUND FRACTIONS 

The systematic uncertainties for the background frac- 
tions are discussed in Ref. [11], and we only summarize 
the values used in this analysis. The systematic uncer- 
tainty on the fraction Jk is set to 9% [11]. The systematic 
uncertainty on the ratio Rk, as indicated in Sec. V, is 
set to half of the uncertainty on ARk given in Eq. (24). 
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FIG. 10: (color online). The ratio Rk obtained using (a) K*'^ 
production, (b) Kg production, and (c) combination of these 
two channels as a function of the kaon transverse momentum. 
The horizontal dashed lines show the mean values. 
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FIG. 12: (color online). The values of (a) Fk, (b) and (c) 
Fp in the like-sign dimuon sample as a function of the kaon, 
pion and proton pt, respectively. The horizontal dashed lines 
show the mean values. 
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FIG. 11: (color online). The difference Rk{K*°) ~ Rk{K%) 
as a function of kaon transverse momentum. The horizontal 
dashed line shows the mean value. 



The systematic uncertainties on the ratios of multiplici- 
ties ni^jriK and n^jriK in pp interactions are set to 4% 
[17]. These multiplicities are required to compute the 
quantities /tt, /p. The ratios Nt^/Nk and Np/Nx, re- 
quired to compute the quantities F^^ and Fp [11] are as- 
signed an additional 4% systematic uncertainty. The val- 
ues of these uncertainties are discussed in Ref. [11]. 



VII. MEASUREMENT OF fs, Fss 

We determine the fraction fs of S muons in the in- 
clusive muon sample and the fraction Fss of events with 
two S muons in the like-sign dimuon sample following 
the procedure described in Ref. [11]. We use the follow- 
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TABLE VII: Asymmetries aK, Itt, and for different pr 
bins. Tfie bottom row shows the mean asymmetries averaged 
over the inclusive muon sample. Only the statistical uncer- 
tainties are given. 



ing value from simulation 

Fll 
FsL + Fll ' 



0.264 ± 0.024, 



(25) 



and obtain 

fs 
Fhkg 
Fll 

FsL 

Fss 



0.536 ± 0.017 (stat) ± 0.043 (syst), 
0.389 ± 0.019 (stat) ± 0.038 (syst), 
0.082 ± 0.005 (stat) ± 0.010 (syst), 

-Fbkg — 2i^LL, 

0.692 ± 0.015 (stat) ± 0.030 (syst). (26) 



The difference between these values and that in Ref. [11] 
are due to the increased statistics and the changes in the 
muon selection and in the analysis procedure. 



VIII. MEASUREMENT OF aK, a„, Up, S 

We measure all detector related asymmetries using the 
methods presented in Ref. [11]. Muons from decays of 
charged kaons and pions and from incomplete absorp- 
tion of hadrons that penetrate the calorimeter and reach 
the muon detectors ("punch-through"), as well as false 
matches of central tracks to segments reconstructed in 
the outer muon detector, are considered as detector back- 
grounds. We use data to measure the fraction of each 
source of background in both the dimuon and inclusive 
muon samples, and the corresponding asymmetries. Data 
are also used to determine the intrinsic charge-detection 
asymmetry of the DO detector. Since the interaction 
length of the meson is greater than that of the K~ 
meson [14], kaons provide a positive contribution to the 
asymmetries A and a. The asymmetries for other back- 
ground sources (pions, protons and falsely reconstructed 
tracks) are at least a factor of ten smaller. 

The results for the asymmetries uk, a^, and Op in dif- 
ferent muon pt bins are shown in Fig. 13 and Table VII. 
The asymmetries and Op are poorly measured in the 
first and last bins due to the small number of events. The 
content of bins 1 and 2, and bins 5 and 6 are therefore 
combined. 



0.08 
0.06 



Bin 


UK X 10^ 


X 10^ 


Up X 10^ 


0.04 


1 
2 


+3.26 ± 1.67 
+4.18 + 0.20 


-0.14 ±0.15 


-6.2 ±6.9 


0.02 


3 


+5.00 ± 0.13 


-0.08 ±0.12 


±4.9 ±5.6 





4 


+5.18 ± 0.22 


±0.25 ± 0.23 


-1.2 ± 12.8 


5 
6 


-h5.44 ± 0.34 
+4.52 ± 0.57 


±0.63 ± 0.40 


-6.8 ±9.6 




All 


+4.88 ± 0.09 


-0.03 ±0.08 


-0.8 ±3.8 
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Pt(K^|x) [GeV] 
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FIG. 13: (color online). The asymmetries (a) aK, (b) a-^, 
and (c) Hp as a function of the kaon, pion and proton px, 
respectively. 



The small residual reconstruction asymmetry Si is mea- 
sured using a sample of J/ip -+ ^J.^^J'~ decays recon- 
structed from two central detector tracks, with at least 
one matching a track segment in the muon detector. The 
values of Si obtained as a function of muon pT are given 
in Table VIII and are shown in Fig. 14. The weighted 
averages for the residual muon asymmetry in the inclu- 
sive muon and the like-sign dimuon samples, calculated 
using weights given by the fraction of muons in each px 
interval {FJ^^) in the inclusive muon (dimuon) sample, 
are given by 



6: 

A ; 



-0.088 ± 0.023)%, 



Eti ^M-^' = (-0.132 ± 0.019)%, 



(27) 
(28) 



where only the statistical uncertainties are given. The 
correlations among different Si are taken into account in 
the uncertainties in Eqs. (27) and (28). 
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TABLE VIII: Muon reconstruction asymmetry Si for different 
muon pt bins. Only the statistical uncertainties are given. 



Bin 



5, X 10^ 



-0.509 ± 
-0.205 ± 
-0.053 ± 
-0.124 ± 
+0.050 ± 
+0.034 ± 



0.106 
0.040 
0.048 
0.075 
0.099 
0.189 



^ 0.005 



0.005 ^ 




15 20 25 
Pt(|x) [GeV] 



FIG. 14: (color online). Muon reconstruction asymmetry as 
a function of muon pT- 



IX. CORRECTIONS FOR BACKGROUND 
ASYMMETRIES 

The corrections for the background and detector con- 
tributions to the measured raw asymmetries a and A are 
obtained combining the resuhs from Tables I, III, VI, and 
VII, and summarized in Tables IX and X. The values in 
the bottom row of these tables are computed by averag- 
ing the corresponding quantities with weights given by 
the fraction of muons in each interval (i^^) in the 
inclusive muon (dimuon) sample, see Eqs. (7) and (8). 
We use the mean values for /^r, F^r, /p, Fp, and Op 
in bins 1 and 2, and in bins 5 and 6, as the number of 
events for those bins are not sufficient to perform sepa- 
rate measurements. 



X. COEFFICIENTS Cb AND Cb 

The dilution coefhcients and Cb in Eq. (12) are ob- 
tained through simulations using the method described 
in Ref. [11]. Both coefficients depend on the value of the 
mean mixing probability, xo- We use the value obtained 
at LEP as averaged by HFAG [3] for this measurement 



Xo(HFAG) = 0.1259 ± 0.0042. 



(29) 



To measure the weights for the different processes pro- 
ducing S muons, we correct the momentum distribution 
of generated h hadrons to match that in the data used in 



TABLE IX: Corrections due to background asymmetries 
/kcik, fwa-w, and fpUp for different pr bins. The bottom row 
shows the weighted averages obtained using weights given by 
the fraction of muons in a given pr interval, in the inclu- 
sive muon sample. Only statistical uncertainties are given. 



Bin 


fxaK X 10^ 


f^a^y X 10^ 


fpap X 10^ 


1 
2 


+0.305 ± 0.220 
+0.624 ± 0.052 


-0.052 ± 0.054 


-0.034 ± 0.041 


3 


+0.832 ± 0.030 


-0.025 + 0.037 


+0.005 + 0.016 


4 


+0.912 ± 0.046 


+0.068 ± 0.065 


-0.008 + 0.081 


5 
6 


+0.785 ± 0.054 
+0.577 ± 0.086 


+0.121 + 0.079 


-0.043 ± 0.077 


All 


+0.776 ± 0.021 


+0.007 + 0.027 


-0.014 + 0.022 


TABLE X: Corrections 


due to background asymmetries 


FxaK 


, fVa,r and FpUp for different pr bins. 


The bottom row 


shows the weighted avera 


ges obtained using 


weights given by 


the fraction of muons in 


a given pT interval, F^, in the like- 


sign dimuon sample. Only statistical uncertainties are given. 


Bin 


FKax X 10^ 


F^a^ X 10^ 


FpUp X 10^ 


1 
2 


+0.300 ± 0.222 
+0.581 ± 0.060 


-0.044 ± 0.046 


-0.029 + 0.035 


3 


+0.719 ± 0.042 


-0.020 ± 0.029 


+0.004 + 0.012 


4 


+0.739 ± 0.050 


+0.050 ± 0.047 


-0.005 + 0.059 


5 
6 


+0.638 ± 0.054 
+0.655 ± 0.112 


+0.094 ± 0.062 


-0.033 ± 0.060 


All 


+0.633 ± 0.031 


-0.002 ± 0.023 


-0.016 + 0.019 



this analysis. The determined weights [17] are given in 
Table XL 

The uncertainty on the weights for the different pro- 
cesses contains contributions from the uncertainty in the 
momentum of the generated b hadrons and from the un- 
certainties in 5-hadron branching fractions. The differ- 
ence in the weights with and without the momentum 
correction contributes to the assigned uncertainties. Ad- 
ditional contributions to the uncertainties on the weights 
derive from the uncertainties on the inclusive branching 



TABLE XI: Heavy-quark decays contributing to the inclusive 
muon and like-sign dimuon samples [17]. The abbreviation 
"non-osc" stands for "non-oscillating," and "osc" for "oscil- 
lating." All weights are computed using MC simulation. 

Process 



Tit 

T2 



fj," X (non-osc) 
h ^ jjT X (osc) 
c ^ ^l+X 

c — >■ ^'^ X (non-osc) 
6 — > c — > 11^ X (osc) 
ccg with c — > fi^ X or c — > /i' 

p", 0(1020), J/i/j, V' ^ At'^M" 
hhcc with c — >■ /i^X or c — >■ X 
cc with c — >■ /i^X or c — > /i^X 



b 
b 

b^ 
b 
b 

b^ 



Wei 


jht 


Wl = 


E 1. 


Wla 


= (1 - xo)wi 


Wlb 


= xo^i^i 


W2 = 


= 0.096 + 0.012 


W2a 


= (1 - X0)UI2 


W2b 




X W3 = 


= 0.064 ± 0.006 


W4 = 


= 0.021 + 0.002 


W5 = 


= 0.013 + 0.002 


We - 


= 0.675 + 0.101 
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TABLE XII: Contribution of different background sources to 
the observed asymmetry in the inclusive muon and like-sign 
dimuon samples. Only statistical uncertainties are given. 



Source 



inclusive muon like-sign dimuon 



{fKttK or FkAk) X lO'' 
iUa^ or F^A^) x 10^ 
{fpttp or FpAp) X 10^ 
[(1 - /bkg)5 or 
(2 - Fbkg)A] X 10^ 
(abkg or ^bkg) X 10^ 



+0.776 ±0.021 
±0.007 ± 0.027 
-0.014 ± 0.022 

-0.047 ± 0.012 



-0.633 ±0.031 
-0.002 ± 0.023 
-0.016 ± 0.019 

-0.212 ± 0.030 



±0.722 ± 0.042 ±0.402 ± 0.053 



(a or A) X 10^ 
[(a — ttbkg) or 
{A - Abkg)] X 10^ 



-0.688 ± 0.002 ±0.126 ±0.041 
-0.034 ± 0.042 -0.276 ± 0.067 



fractions B fiX, B ^ cX and B ^ cX [14]. We as- 
sign an additional uncertainty of 15% to the weights 
and uig for uncertainties on the cross sections for cc and 
bbcc production. 

The resulting ct, and Cb coefficients are found to be 



Q, = ±0.061 + 0.007, 
C'b = ±0.474 + 0.032. 

XI. ASYMMETRY A^i 



(30) 
(31) 



The results obtained in Sees. IV-X are used to measure 
the asymmetry following the procedure of Ref. [11]. 
Using 2.041 x 10*^ muons in the inclusive muon sample 
and 6.019 X 10® events in the like-sign dimuon sample we 
obtain the following values for the uncorrected asymme- 
tries a and A: 



a = (+0.688 + 0.002)%, 
A = (+0.126 + 0.041)%. 



(32) 
(33) 



The difference between these values and those in Ref. [11] 
are due to increased statistics and the changes in the 
muon selection. The contributions from different back- 
ground sources to the observed asymmetries a and A are 
summarized in Table XII. 

The asymmetry A^^, extracted from the asymmetry a 
of the inclusive muon sample using Eqs. (7) and (30), is 



(-1.04 + 1.30 (stat) + 2.31 (syst))%. (34) 



The contributions to the uncertainty are given in Ta- 
ble XIII. Figure 15(a) shows a comparison of the asym- 
metry a and the background asymmetry, Obkg = fs^ + 
fxciK ± /ttOtt ± fpO-p, as a function of muon px- There 
is excellent agreement between these two quantities, with 
X^/d.o.f. = 0.8/6 for their difference. Figure 15(b) shows 
the value of /sag ~ a — flbkg, which is consistent with 
zero. The values a and abkg are given in Table XIV. This 
result agrees with the expectation that the value of the 
asymmetry a is determined mainly by the background, 
as the contribution from is strongly suppressed by 



TABLE XIII: Sources of uncertainty on A^i from Eqs. (34), 
(35), and (36). The first nine rows contain statistical uncer- 
tainties, while the next four rows reflect contributions from 
systematic uncertainties. 



Source 



HA'^i) X 10^ 



5{At,) X 10^ 



S{Asi) X 10^ 





-Cjq. 


ii;q. \^oO } 


iijq. yoo ) 


A 7 \ 

A OY a (stat) 


n nan 


n 1 01 
U.izi 


O 1 Q9 




n 479 




n n9s 




M / A 


u.uoy 


U.UDO 




n 1 SI 


u.uzo 


u.uuo 




0.323 


0.026 


0.002 


Ak 


0.458 


0.052 


0.037 




0.802 


0.067 


0.030 


Ap 


0.584 


0.050 


0.020 


5 or A 


0.377 


0.087 


0.067 


Ik (syst) 


2.310 


0.204 


0.007 


Rk (syst) 


N/A 


0.068 


0.072 


n, K, p multiplicity 


0.067 


0.019 


0.017 


Cb or Cb 


0.121 


0.052 


0.056 


Total statistical 


1.304 


0.202 


0.172 


Total systematic 


2.313 


0.222 


0.093 


Total 


2.(i5() 


0.:!00 


0.196 



TABLE XIV: The measured asymmetry a and the expected 
background asymmetry abkg in the inclusive muon sample for 
different pr bins. For the background asymmetry, the first 
uncertainty is statistical, the second is systematic. 



bin 



a X 10^ 



Obkg X 10^ 



-0.071 ± 

±0.503 ± 
±0.712 ± 
±0.841 ± 
±0.812 ± 
+0.702 ± 



0.025 
0.005 
0.003 
0.005 
0.007 
0.010 



-0.055 ± 
+0.438 ± 
+0.785 ± 
±0.910 ± 
±0.897 ± 
+0.680 ± 



0.240 ± 0.664 
0.089 + 0.117 
0.056 + 0.083 
0.124 + 0.105 
0.139 + 0.101 
0.189 ± 0.059 



the factor of Cb = 0.061 + 0.007. The consistency of 
with zero in Eq. (34) and the good description of 
the charge asymmetry a for different values of muon pt 
shown in Fig. 15 constitute important tests of the validity 
of the background model and of the method of analysis 
discussed in this Article. 

The second measurement of the asymmetry A^i, ob- 
tained from the uncorrected asymmetry A of the like-sign 
dimuon sample using Eqs. (8), (30) and (31), is 



(-0.808 + 0.202 (stat) + 0.222 (syst))%, (35) 



where we take into account that both as and As in 
Eq. (8) are proportional to A^j, and that FssCb+FsLCb = 
0.342 + 0.028. The contributions to the uncertainty of A^^ 
for this measurement are also listed in Table XIII. 

The measurement of the asymmetry A'^^ using the lin- 
ear combination given in Eq. (13) is performed following 
the procedure described in Ref. [11]. We select the value 
of the parameter a that minimizes the total uncertainty 
on the measurement. Appendix A gives more de- 
tails on this method of combination. All uncertainties in 
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FIG. 15: (color online), (a) The asymmetry flbkg (points with 
error bars representing the total uncertainties) , expected from 
the measurements of the fractions and asymmetries for back- 
ground processes, is compared to the measured asymmetry a 
for the inclusive muon sample (shown as a histogram, since 
the statistical uncertainties are negligible). The asymmetry 
from CP violation is negligible compared to the background 
in the inclusive muon sample, (b) The difference a — Obkg. 
The horizontal dashed line shows the mean value. 



Table XIII, except the statistical uncertainties on a. A, 
and Rk, are treated as fully correlated. This leads to 
a = 0.89, and the corresponding value of the asymmetry 
is 



A'' 



4^ = (-0.787 ± 0.172 (stat) ± 0.093 (syst))%. (36) 

This value is used as the final result for A^j. It differs by 
3.9 standard deviations from the standard model predic- 
tion of given in Eq. (5). The different contributions 
to the total uncertainty on A^j in Eq. (36) are listed in 
Table XIII. 

The measured value of A^j places a constraint on the 
charge asymmetries a^j and a|j. The asymmetry A^j is a 
linear combination of the semi-leptonic charge asymme- 
tries from and B° meson decays [2]. The coefficients 
Cd and Cs in Eq. (2) depend on the mean mixing prob- 
ability and the production rate of -B" and B'^ mesons. 
We use the latest measurements of these quantities from 
LEP as averaged by HFAG [3] 

Xo(HFAG) = 0.1259 ± 0.0042, (37) 
/d(HFAG) = 0.403 ±0.009, (38) 
/s(HFAG) = 0.103± 0.009, (39) 



FIG. 16: (color online). Comparison of in data with the 
SM prediction for afi and afi. Also shown are the measure- 
ments of flgi [3] and a|i [19]. The error bands represent the 
±1 standard deviation uncertainties on each individual mea- 
surement. The 95% C.L. band is also given for this mea- 
surement. 



and find the values given in Eq. (4). 

Figure 16 presents the measurement in the (afi,agi) 
plane together with the existing direct measurements of 
a^j from the B factories [3] and the independent DO mea- 
surement of a|j in B'^ — fJ-DgX decays [19]. All measure- 
ments are consistent. 

The quantity Aj-cs defined as 



Ar 



{A 



bkg 



aabkg) 



(40) 



is the residual charge asymmetry of like-sign dimuon 
events after subtracting all background contributions 
from the raw charge asymmetry. This quantity does not 
depend on the interpretation in terms of the charge asym- 
metry of semi-leptonic decays of B mesons. We obtain 

Arcs = (-0.246 ± 0.052 (stat) ± 0.021 (syst))%, (41) 

The measured value of A^cs differs by 4.2 standard devi- 
ations from the standard model prediction 



ylros(SM) = (-0.009 ± 0.002)%. 



XII. CONSISTENCY CHECKS 



(42) 



To study the stability of the result, we repeat this 
measurement with modified selections, and with subsets 
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of the available data. The only difference compared to 
Ref. [11] is Test D, where we applied a stronger criterion 
on the muon IP, following the suggestion of Ref. [18] . In 
all tests the modified selections were applied to all muons. 
For completeness, we give the full list of tests performed: 

• Test Al: Using only the part of the data sample 
corresponding to the first 2.8 fb~^. 

• Test A2: Using only the part of the data sample 
corresponding to the previous measurement with 

6.1 fb-i [11]. 

• Test A3: Using only the part of the data sample 
corresponding to the last 2.9 fb~^. 

• Test B: In addition to the reference muon selections 
[11], we require at least three hits in the muon wire 
chambers (layers B or C), and lower the require- 
ment for the fit to a track segment reconstructed 
in the muon detector. 

• Test C: Since background muons are mainly pro- 
duced by decays of kaons and pions, their track 
parameters measured in the central tracker and by 
the muon system can differ. The background frac- 
tion therefore depends strongly on the of the 
difference between these two measurements. The 
requirement on this is changed from 12 to 4. 

• Test D: The maximum value of the IP is changed 
from 0.3 to 0.012 cm. This test is also sensitive to 
possible contamination from cosmic-ray muons. 

• Test E: Using low-luminosity data with fewer than 
three interaction vertices. 

• Test F: Using events corresponding to only two of 
four possible configurations for the magnets, with 
identical solenoid and toroid polarities. 

• Test G: Changing the minimum requirement on the 
invariant mass of the two muons from 2.8 GeV to 
12 GeV. 

• Test H: Using the same muon -pT requirement, pt > 

4.2 GeV, for the full detector acceptance. 

• Test I: Requiring the muon to be < 7.0 GeV. 

• Test J: Requiring the azimuthal angle <\) of the muon 
track to be in the range < < 4 or 5.7 < (f) < 27r. 
This selection excludes muons with reduced muon 
identification efficiency in the region of the support 
structure of the detector. 

• Test K: Requiring the muon rj to be in the range 
|?7| < 1.6. This test is sensitive to possible contam- 
ination from muons associated with beam halos. 

• Test L: Requiring the muon r} to have |r?] < 1.2 or 
1.6 < |r/| < 2.2. 



• Test M: Requiring the muon rj to be in the range 
|77| < 0.7 or 1.2 < \ri\ < 2.2. 

• Test N: Requiring the muon t] to be in the range 
0.7 < \r]\ < 2.2. 

• Test 0: Using like-sign dimuon events that pass at 
least one single muon trigger, while ignoring the 
requirement for a dimuon trigger. 

• Test P: Using like-sign dimuon events passing both 
single muon and dimuon triggers. 

A summary of these studies is presented in Tables XV 
and XVI. The last row, denoted as "Significance" , gives 
the absolute value of the difference between the reference 
result (column Ref) and each modification, divided by 
its uncertainty, taking into account the overlap in events 
between the reference and test samples. Both statisti- 
cal and systematic uncertainties are used in the calcula- 
tion of the significance of the difference. The x^ of these 
tests defined as the sum of the square of all significances 
is x^ = 17.1 for 18 tests. These tests demonstrate the 
stability of the measured asymmetry A^j, and provide a 
confirmation of the validity of the method. 

We also compare the dependence on the muon pseudo- 
rapidity T]{ii) of the observed and expected charge asym- 
metry in the inclusive muon sample. We repeat the anal- 
ysis procedure, but measure all background contributions 
as a fimction of |?7(/i)|. The result of this comparison 
is shown in Fig. 17. The dependence on |r?(/Lt)| is cor- 
rectly described by the background asymmetry. There 
is good agreement between these two quantities, with a 
X^/d.o.f. = 2.8/4. This is consistent with our expecta- 
tion that the contribution of A^j in the inclusive muon 
charge asymmetry is overwhelmed by background. 

Figure 18 shows the observed and expected uncor- 
rected like-sign dimuon charge asymmetry as a function 
of the dimuon invariant mass. The expected asymmetry 
is computed using Eq. (8) and the measured parameters 
of sample composition and asymmetries. As in Ref. [11], 
we compare the expected uncorrected asymmetry using 
two assumptions for A^j. In Fig. 18(a) the observed asym- 
metry is compared to the expectation for the SM value 
of A^i(SM) = -0.028%, while Fig. 18(b) shows the ex- 
pected asymmetry for A^^ = —0.787%. Large discrepan- 
cies between the observed and expected asymmetries can 
be observed for = Agi(SM), while good agreement 
is obtained for the measured A^^ value corresponding to 
Eq. (36). The observed asymmetry changes as a function 
of dimuon invariant mass, and the expected asymmetry 
tracks this effect when A^^ = —0.787%. This dependence 
of the asymmetry on invariant mass of the muon pair is a 
function of the production mechanism of the particles in- 
volved and of their decays. The agreement between the 
observed and expected asymmetries indicates that the 
physics leading to the observed asymmetry can be de- 
scribed by contributions from the background and from 
decays of b hadrons. 
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TABLE XV: Measured asymmetry A^i for the reference selection (column Ref) and for samples used in Tests A - H. 







Ref 


Al 


A2 


A3 


B 


C 


D 


E 


F 


G 


H 


Niii/i) X 10" 
a X 10^ 


6 


6.019 


1.932 


3.991 


2.028 


4.466 


3.280 


2.857 


3.128 


3.012 


2.583 


2.220 




+0.688 


+0.703 


+0.680 


+0.702 


+0.548 


+0.325 


+0.835 


+0.682 


+0.727 


+0.688 


+0.751 


A X 10^ 




+0.126 


+0.061 


+0.062 


+0.259 


-0.149 


-0.361 


+0.555 


+0.136 


+0.137 


+0.450 


+0.344 


a 




0.894 


0.760 


0.851 


0.813 


0.891 


0.631 


1.271 


0.831 


0.940 


0.939 


0.807 


[(2 - Fbkg)A 


- afsS] X 10^ 


-0.170 


-0.193 


-0.178 


-0.157 


-0.270 


-0.370 


-0.133 


-0.206 


-0.152 


-0.114 


-0.049 


fs 




0.536 


0.583 


0.557 


0.516 


0.509 


0.560 


0.472 


0.534 


0.493 


0.537 


0.536 


Jbkg 




0.389 


0.336 


0.365 


0.384 


0.405 


0.338 


0.627 


0.374 


0.407 


0.436 


0.325 


A^i X 10^ 




-0.787 


-0.803 


-0.891 


-0.600 


-0.906 


-0.708 


-1.138 


-0.584 


-0.986 


-0.379 


-0.654 


a(^Ji) X 10" 


(stat) 


0.172 


0.278 


0.204 


0.335 


0.207 


0.220 


0.365 


0.224 


0.302 


0.263 


0.254 


a(^'i) X 10^ 


(syst) 


0.093 


0.125 


0.128 


0.188 


0.107 


0.104 


0.323 


0.108 


0.135 


0.209 


0.103 


Significance 






0.007 


0.742 


0.567 


1.029 


0.525 


1.022 


1.236 


0.960 


1.537 


1.120 



TABLE XVI: Measured asymmetry A^i for the reference selection (column Ref) and for samples used in Tests I - P. 







Ref 


1 


J 


K 


L 


M 


N 





P 


7V(/i/i) X 10" 


t> 


6.019 


4.428 


3.504 


2.928 


2.741 


4.259 


3.709 


2.724 


2.440 


ax 10^ 




+0.688 


+0.672 


+0.691 


+0.711 


+0.761 


+0.501 


+0.802 


+0.688 


+0.688 


A X 10^ 




+0.126 


+0.250 


+0.160 


+0.118 


+0.216 


-0.033 


+0.262 


+0.245 


+0.272 


a 




0.894 


0.908 


0.817 


0.872 


0.825 


0.702 


0.908 


0.941 


0.898 


[(2 - Fbkg)A 


- afsS] X 10" 


-0.170 


-0.209 


-0.187 


-0.221 


-0.214 


-0.187 


-0.150 


-0.126 


-0.122 


fs 




0.536 


0.514 


0.555 


0.556 


0.570 


0.519 


0.514 


0.536 


0.536 


-Fbkg 




0.389 


0.414 


0.352 


0.363 


0.333 


0.402 


0.428 


0.408 


0.395 


A^i X 10" 




-0.787 


-0.925 


-0.569 


-0.847 


-0.430 


-0.761 


-0.774 


-0.809 


-0.689 


a(^^) X 10^ 


(stat) 


0.172 


0.204 


0.202 


0.224 


0.260 


0.207 


0.221 


0.247 


0.253 


a{Ai,) X 10" 


(syst) 


0.093 


0.115 


0.100 


0.122 


0.117 


0.110 


0.118 


0.129 


0.128 


Significance 






1.245 


1.672 


0.377 


1.678 


0.441 


0.186 


0.120 


0.497 



We also measure the mean mixing probability using the 
ratio of like-sign and opposite-sign dimuon events. The 
background contribution in both samples is obtained us- 
ing the method presented in this Article. The measured 
mean mixing probability is found to be consistent with 
the world average value [3] . 

We conclude that our method of analysis provides a 
consistent description of the dimuon charge asymmetry 
for a wide range of input parameters, as well as for sig- 
nificantly modified selection criteria. 

XIII. COMPARISON WITH SIMULATION 

The measurement of the background fractions is based 
on data, and the input from simulation is limited to the 
ratio of multiplicities ri^/n^ and np/riK in pp interac- 
tions [11]. Nevertheless, it is instructive to compare the 
results obtained in data and in simulation. Such a com- 
parison is shown in Table XVII. The simulation used 
in this analysis is described in Ref. [11]. All quantities 
measured in simulation are obtained using the informa- 
tion on the generated processes. All uncertainties in the 
second and third columns are statistical. The difference 
between the values obtained in data and simulation is 
given in the fourth column and includes the systematic 
uncertainties. The agreement between the measured and 



simulated quantities is satisfactory. The excellent agree- 
ment between the mean values of Rk, which is one of 
the most essential quantities of this measurement and for 
which many systematic uncertainties cancel, is especially 
notable: 

i?K(data) = 0.856 ±0.020 (stat) ±0.026 (syst), 
Rk{MC) = 0.901 ± 0.086 (MC stat). (43) 

This comparison provides support for the validity of the 
presented measurement. 

XIV. DEPENDENCE OF ASYMMETRY A^i ON 
MUON IMPACT PARAMETER 

The asymmetry A^^ is produced by muons from direct 
semi-leptonic decays of b quarks. A distinctive feature of 
these muons is the large impact parameter of their tra- 
jectories with respect to the primary vertex [12, 18]. The 
simulation shows that the dominant source of background 
from L muons corresponds to charged hadrons produced 
in the primary interactions that then decay to muons, 
and the tracks of such muons have small impact parame- 
ters if the decay is outside the tracking volume. Figure 19 
shows the muon IP distribution in data and in simula- 
tion. The shaded histogram shows the contribution from 
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FIG. 17: (color online), (a) The asymmetry flbkg (points with 
error bars representing the total uncertainties), as expected 
from the measurements of the fractions and asymmetries of 
the background processes, is compared to the measured asym- 
metry a of the inclusive muon sample (shown as a histogram, 
since the statistical uncertainties are negligible) as a func- 
tion of the absolute value of muon pseudorapidity \rj{ii)\. The 
asymmetry from CP violation is negligible compared to the 
background in the inclusive muon sample, (b) The difference 
a — fflbkg. The horizontal dashed line shows the mean value. 



TABLE XVII: Comparison of background fractions measured 
in data and in simulation. Only the statistical uncertainties 
are given in the second and third column. The difference 
between data and simulation is given in the fourth column 
and includes both statistical and systematic uncertainties. 



Quantity 


Data 


Simulation 


Difference 


fK X 10" 


15.96 ± 0.24 


14.31 ± 0.06 


-f 1.65 ± 2.55 


X 10^ 


30.01 ± 1.60 


29.82 ± 0.09 


+0.19 ± 5.15 




0.38 ± 0.17 


1.07 ± 0.02 


-0.69 ± 0.60 


Fk X 10^ 


13.78 ± 0.42 


12.89 ± 1.32 


-hO.89 ± 2.26 


F^ X 10^ 


24.81 ± 1.38 


25.88 ± 1.86 


-1.07 ± 4.36 


Fp X 10^ 


0.35 ± 0.14 


1.29 ± 0.39 


-0.94 ± 0.72 


fs X 10" 


53.65 ± 1.74 


54.79 ± 0.14 


-1.14 ± 7.11 


Fbkg X lO"* 


38.94 ± 1.89 


40.01 ± 2.31 


-1.07 ± 6.21 


Rk X 10^ 


85.62 ± 1.98 


90.08 ± 8.60 


-4.46 ± 9.74 



L muons in simulation, which decreases significantly for 
increasing values of the muon IP. The background can 
therefore be significantly suppressed by selecting muons 
with large impact parameter. 

To verify the origin of the observed charge asymmetry, 
we perform several complementary measurements. We 



FIG. 18: (color online). The observed and expected like-sign 
dimuon charge asymmetries in bins of dimuon invariant mass. 
The expected asymmetry is shown for (a) = ^^[(SM) and 
(b) Al^ = -0.787%. 
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FIG. 19: (color online). The muon IP distribution in the 
inclusive muon sample (bullets) . The solid line represents the 
muon IP distribution in simulation. The shaded histogram is 
the contribution from L muons in simulation. 



require the muon IP to be larger or smaller than 120 fini. 
For events in the like-sign dimuon sample, we require that 
both muons satisfy these conditions. These measure- 
ments are denoted as /P>i2o and /P<i2o, respectively. 
The selected threshold of 120 /im can be compared with 
the spread in the crossing point of the colliding beams in 
the Tevatron collider, and with the precision of pp vertex 
reconstruction, which are about 30 /im and 15 fxm, re- 
spectively, in the plane perpendicular to the beam axis. 
The chosen value of 120 fim gives the minimal uncer- 
tainty on Ogj and defined in Eq. (2). 
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TABLE XVIII: Input quantities for the measurement of A^i 
using rnuons with IP above and below 120 //m. Only statis- 
tical uncertainties are given. 



Quantity 


IP > 120 urn 


IP < 120 ^im 


fK X 10^ 


5.19 ± 0.37 


17.64 ± 0.27 


U X 10^ 


5.65 ± 0.40 


34.72 ± 1.86 


fp X 10^ 


0.05 ± 0.03 


0.45 ± 0.20 


Fk X 10^ 


4.48 ± 4.05 


21.49 ± 0.62 


X 10^ 


4.43 ± 3.95 


40.47 ± 2.26 


Fp X 10^ 


0.03 ± 0.05 


0.59 ± 0.23 


fs X 10" 


89.11 ± 0.88 


47.18 ± 2.03 


Fbkg X 10^ 


8.94 ± 8.26 


62.56 ± 3.07 


Fss X 10^ 


91.79 ± 7.65 


53.66 ± 2.68 


W A ±U 






abkg X 10^ 


+0.027 ± 0.023 


+0.864 + 0.049 


A X 10^ 


-0.529 ± 0.120 


+0.555 ± 0.060 


.4bkg X 10^ 


-0.127 ± 0.093 


+0.829 ± 0.077 




0.70 ± 0.05 


0.95 ± 0.02 


Ck 


0.39 ± 0.06 


0.98 ± 0.01 


Fll/{Fll + Fsl) 


0.089 ± 0.062 


0.350 ± 0.029 


Cb 


0.109 ± 0.011 


0.038 ± 0.007 


Cb 


0.526 ± 0.037 


0.413 ± 0.032 



In total, 0.356 x 10'^ muons in the inclusive muon sam- 
ple and 0.714 x 10^ events in the like-sign dimuon sam- 
ple are selected for the /P>i2o measurement. Events 
are subject to the same analysis as for the entire sam- 
ple, except that the ratio RK(Kg) is not used because of 
insufficient Kg — > tt+tt" decays in the dimuon sample. 
Background asymmetries should not depend on the muon 
IP, and wc verified that the difference in kaon asymme- 
try for the whole sample and the /P>i2o events agree: 
aif(/P>i2o) - aif(all) = (-1.6 ± 1.5)%. We therefore 
use the values given in Tables VII and VIII. All other 
measured quantities are given in Table XVIII. The back- 
ground fractions are strongly suppressed in the /P>i20 
sample, and their influence on the measurement of is 
significantly smaller. Using these values, we obtain for 
the inclusive muon sample 

A^i(7P>i2o) = (-0.422 ± 0.240 (stat) ± 0.121 (syst))%, 

and for the like-sign dimuon sample 

-4si(/P>i2o) = (-0.818 ± 0.342 (stat) ± 0.067 (syst))%. 

(45) 

We obtain the final value of ^g[(/P>i2o) using the lin- 
ear combination of Eq. (13), and select the value of a to 
minimize the total uncertainty on A'^^, which corresponds 
to a — —9.29. The combination for a negative value of a 
is equivalent to the weighted average of Eqs. (44) and (45) 
taking into account the correlation of uncertainties (see 
Appendix A for more details). The corresponding asym- 
metry is found to be 

^^i(/P>i2o) = (-0.579 ± 0.210 (stat) ± 0.094 (syst))%. 

(46) 



TABLE XIX: Sources of uncertainty on A'^i{IPyi2o) in 
Eqs. (44), (45), and (46). The first nine rows contain statis- 
tical uncertainties, and the next four rows contain systematic 
uncertainties. 



Source 


(5(4,^) X 10^ 


S(Ali) X 10^ S{A'',i) X 10- 




Eq. (44) 


En (45) 


En (46) 


^ or a (stat) 


0.055 


0.244 


0.093 


Jk (stat) 


0.048 


0.031 


0.058 


Rk (stat) 


N/A 


0.244 


0.074 


Pin fi)/P{K fi) 


0.007 


0.004 


0.006 


P{p n)/P{K fi) 


0.012 


0.004 


0.010 


Ak 


0.023 


0.012 


0.017 




0.037 


0.009 


0.026 


Ap 


0.025 


0.007 


0.019 


(5 or A 


0.210 


0.075 


0.157 


/if (syst) 


0.112 


0.027 


0.083 


Rk (syst) 


N/A 


0.014 


0.007 


n, K, p multiplicity 


0.016 


0.016 


0.016 


Cb or Cb 


0.043 


0.057 


0.041 


Total statistical 


0.240 


0.342 


0.210 


Total systematic 


0.121 


0.067 


0.094 


Total 


0.269 


0.348 


0.230 



The contributions to the uncertainties in Eqs. (44 - 46) 
arc given in Table XIX. 

From the known frequencies of oscillations, AMq/2n 
{q = d,s), the period of oscillation for the i?'^ meson is 
many times longer than its lifetime so that the mixing 
probability of P" mesons effectively increases with long 
decay lengths and large impact parameters. The P° me- 
son oscillates a number of times within its lifetime so that 
it is "fully mixed" for any appreciable impact parameter 
requirement. As a result, the fraction of P" mesons that 
have oscillated into the other fiavor is increased in the 
sample with large muon impact parameter. This behav- 
ior is demonstrated in Fig. 20, which shows the normal- 
ized IP distributions for muons produced in oscillating 
decays of P*^ and P° mesons in simulation. The contri- 
bution of the ttgj asymmetry in is therefore enhanced 
in the sample with a large muon IP. From simulation, the 
mixing probability of P° meson in the /P>i20 sample is 
determined to be 

Xd(/P>i20,MC) = 0.342+0.004, (47) 

with the uncertainty limited by the number of simulated 
events. This value can be compared to the input to the 
simulation for the P° mixing probability integrated over 
time, Xd = 0.1864 + 0.0022 [3]. The coefficients Cd and 
Cs in Eq. (2) for the /P>i2o selection become 

Cd{IP>i2o) = 0.728 ±0.018, 

Cs{IP>i2o) = 0.272 + 0.018. (48) 

The value of Ag[(/P>i2o) should therefore be reduced rel- 
ative to the value for the full dimuon sample, if the con- 
tribution of afj dominates the asymmetry A^^. 

The measurement of /P<i2o is performed using 1.687 x 
10^ muons in the inclusive muon sample and 2.857 x 10^ 
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FIG. 20: (color online). The normalized IP distribution for 
muons produced in oscillating decays of mesons (points 
with error bars) and mesons (solid histogram) in simula- 
tion. 



TABLE XX: Sources of uncertainty on y4'i(/P<i2o) in 
Eqs. (49), (50), and (51). The first nine rows contain sta- 
tistical uncertainties, the next four rows contain systematic 
uncertainties. 



Source 



5{A'^i) X 10^ S{A^^i) X 10^ 



lO'' 





Eq. (49) 


Eq. (50) 


Eq. (51) 


A or a (stat) 


0.136 


0.233 


0.285 


/k (stat) 


1.059 


0.173 


0.082 


Rk (stat) 


N/A 


0.141 


0.155 


P{7V ^ ^)/P{K ^ /i) 


0.388 


0.060 


0.026 


P{p ^ fi)/P{K ^ /i) 


0.699 


0.064 


0.004 


Ak 


0.986 


0.123 


0.089 


Ati 


1.727 


0.165 


0.075 


Ap 


1.261 


0.123 


0.050 


(5 or A 


0.606 


0.107 


0.071 


fx (syst) 


4.951 


0.508 


0.034 


Rk (syst) 


N/A 


0.286 


0.307 


n, K, p multiplicity 


0.137 


0.034 


0.025 


Cb or Cb 


0.305 


0.087 


0.093 


Total statistical 


2.774 


0.439 


0.366 


Total systematic 


4.962 


0.590 


0.323 


Total 


5.685 


0.735 


0.488 



events in the like-sign dimuon sample. Exactly the same 
procedure is applied as for the main measurement, using 
the background and muon reconstruction asymmetries 
given in Tables VII and VIII. All other quantities are 
given in Table XVIII. The background fractions are sig- 
nificantly increased in the samples with small muon IP, 
thereby increasing the uncertainties related to the back- 
ground description (Table XX). 

Using these values we obtain from the inclusive muon 
sample 

AliiIP<i2a) = (-1.65±2.77 (stat)±4.96 (syst))%, (49) 

and from the like-sign dimuon sample 

^si(^-P<i2o) = (-1.17±0.44 (stat)±0.59 (syst))%. (50) 

The measurement using the linear combination given in 
Eq. (13) is performed with a — +1.27, which minimizes 



the total uncertainty on A^^. The value of A^^ is found to 
be 

-4^i(/^'<i2o) (-1.14±0.37 (stat)±0.32 (syst))%. (51) 

The mean mixing probability Xd in the /P<i2o sample 
obtained in simulation is found to be 



Xd{IP<i2o,MC) = 0.084 ± 0.002, 



(52) 



and the coefhcients Cd and Cg in Eq. (2) for the /P<i2o 
selection are 



Cd{IP<120) 
Cs{IP< 120) 



= 0.397 ±0.022, 
= 0.603 ±0.022. 



(53) 



The measurements with //<i2o and /P>i2o use in- 
dependent data samples, and the dependence of A^^ on 
flgj and is different for the /P<i2o and /P>i2o sam- 
ples. The measurements given in Eqs. (46) and (51) can 
therefore be combined to obtain the values of a^j and 
a|[, taking into account the correlation among different 
sources of uncertainty. All uncertainties in Tables XIX, 
XX, except the statistical uncertainties on a, A, fx, Rk, 
P(7r ^ fi)/P{K n), and P{p n)IP{K are 
treated as fully correlated. The values of a^j and a^j ex- 
tracted are 

afi = (-0.12 ±0.52)%, 

<i = (-1.81 ± 1.06)%. (54) 

The correlation pds between these two quantities is 



Pds 



-0.799. 



(55) 



The uncertainty on a^j and a^j obtained in this study is 
comparable with that obtained from the direct measure- 
ments. Figure 21 presents the results of the IP study in 
the (asi,afi) plane together with the result (36) of the 
measurement using all like-sign dimuon events. The 
ellipses represent the 68% and 95% two-dimensional C.L. 
regions, respectively, of a|j and a^j values obtained from 
the measurements with IP selections. 

We also performed four additional measurements with 
IP thresholds of 50 /xm and 80 /xm. They are denoted 
as IP<5Q, IP>50, IP<80, and IP>so, respectively. The 
input quantities for these measurements are presented in 
Tables XXI and XXII. The A^j values in the inclusive 
and like-sign dimuon samples and their combinations are 
given in Table XXIII. The mean mixing probability Xd 
for all these measurement is obtained through simulation. 
The results are presented in Table XXIV, together with 
the corresponding coefhcients Cd and Cs ■ 

As for the combinations of the /P<i2o and /P>i2o sam- 
ples, the measurements with /P<5o and /P>5o samples, 
as well as with /P<80 and /P>80 samples, can be com- 
bined to determine the values of a^j and (Table XXV). 
The measurements with different IP thresholds are con- 
sistent with each other within two standard deviations 
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FIG. 21: (color online). Measurements of with different 
muon IP selections in the (afi, al.-^ plane. The bands represent 
the ±1 standard deviation uncertainties on each individual 
measurement. The ellipses represent the 68% and 95% two- 
dimensional C.L. regions, respectively, of afi and a^i values 
obtained from the measurements with IP selections. 



taking into account the correlation between the uncer- 
tainties. 

We conclude that the observed dependence of the like- 
sign dimuon charge asymmetry on muon IP is consis- 
tent with the hypothesis that it has its origin from semi- 
leptonic 6-hadron decays. The contributions of a^j and 
flgj to can be determined separately by dividing the 
sample according to the muon IP, although the uncer- 
tainties on the values of a^j and do not allow for the 
definitive conclusion that the deviation of from its 
SM prediction is dominated from the a|[ asymmetry. 



XV. CONCLUSIONS 

We have presented an update to the previous measure- 
ment [11] of the anomalous like-sign dimuon charge asym- 
metry with 9.0 fb"-'^ of integrated luminosity. The 
analysis has improved criteria for muon selection, which 
provide a stronger background suppression and increase 
the size of the like-sign dimuon sample. A more accu- 
rate measurement of the fraction of kaons that produce 
muons in the inclusive muon sample (fx), and an addi- 
tional measurement of the ratio of such yields in like-sign 
dimuon to inclusive muon data {Rk = Fx/ fK) using 
— > 7r''"7r~ decay have been performed. This provides 
better precision of Rk, and an independent estimate of 



TABLE XXI: Input quantities for the measurement of A^i 
using muons with IP above 50 ^m, 80 /xm and 120 /^m, re- 
spectively. Only statistical uncertainties are given. 



Quantity 


IP > 50 ^im 


IP > 80 fj,m 


IP > 120 ^m 


fx X 10" 


6.47 ± 0.18 


5.38 ± 0.24 


5.19 ± 0.37 


X 10^ 


10.42 ± 0.47 


7.24 ± 0.38 


5.65 ± 0.40 


U X 10^ 


0.11 ± 0.05 


0.07 ± 0.03 


0.05 ± 0.03 


Fk X 10^ 


6.31 ± 1.73 


4.79 ± 2.59 


4.48 ± 4.05 


X 10^ 


9.51 ± 2.36 


6.39 ± 2.95 


4.43 ± 3.95 


Fp X 10^ 


0.11 ± 0.06 


0.03 ± 0.04 


0.03 ± 0.05 



fs X 10-^ 
i^bkg X 10" 
Fss X 10^ 



82.99 
15.91 
85.63 



± 0.81 
± 4.38 
± 3.74 



87.32 
11.39 
89.88 



± 0.74 
± 6.10 
± 5.10 



89.11 
8.94 
91.79 



± 0.88 
± 8.26 
± 7.65 



a X 10" 
dbkg X 10^ 
A X 10^ 
^bkg X 10^ 



+0.134 
+0.146 
-0.302 
-0.043 



± 0.004 
± 0.024 
± 0.079 
± 0.071 



+0.035 
+0.068 
-0.386 
-0.139 



± 0.005 
± 0.023 
± 0.094 
± 0.083 



-0.014 
+0.027 
-0.529 
-0.127 



± 0.005 
± 0.023 
± 0.120 
± 0.093 



Ck 



Fll 



(Fll+Fsl) 
Cb 



0.81 
0.66 
0.108 
0.084 
0.496 



± 0.03 
± 0.03 
± 0.038 
± 0.008 
± 0.034 



0.75 
0.52 
0.125 
0.095 
0.510 



± 0.05 
± 0.05 
± 0.060 
± 0.009 
± 0.034 



0.70 
0.39 
0.089 
0.109 
0.526 



± 0.05 
± 0.06 
± 0.062 
± 0.011 
± 0.037 



TABLE XXII: Input quantities for the measurement of 
using muons with IP below 50 fim, 80 fim and 120 fim, re- 
spectively. Only statistical uncertainties are given. 



Quantity 


IP < 50fim 


IP < 80 fj,m 


IP < 120 ^m 


fx X 10" 


19.35 ± 0.33 


18.32 ± 0.30 


17.64 ± 0.27 


.fn X 10^ 


37.58 ± 2.08 


34.34 ± 1.95 


34.72 ± 1.86 


.fp X 10^ 


0.51 ± 0.22 


0.48 ± 0.21 


0.45 ± 0.20 


Fk X 10^ 


28.03 ± 0.95 


23.79 ± 0.74 


21.49 ± 0.62 


X 10^ 


51.72 ± 3.18 


44.26 ± 2.63 


40.47 ± 2.26 


Fp X 10^ 


0.77 ± 0.29 


0.66 ± 0.25 


0.59 ± 0.23 



fs X 10" 
Fbkg X 10" 
Fss X 10^ 



42.56 
81.53 
43.42 



± 2.73 
± 4.30 
± 3.75 



45.40 
70.13 
48.76 



± 2.13 
± 3.52 
± 2.84 



47.18 
62.56 
53.66 



± 2.03 
± 3.07 
± 2.68 



a X 10" 
dbkg X 10^ 
A X 10^ 



^bkg 



X 10' 



+0.953 
+0.997 
+0.715 
+1.243 



± 0.003 
± 0.056 
± 0.083 
± 0.096 



+0.896 
+0.916 
+0.683 
+0.994 



± 0.003 
± 0.052 
± 0.069 
± 0.082 



+0.835 
+0.864 
+0.555 
+0.829 



± 0.002 
± 0.049 
± 0.060 
± 0.077 



Ck 

Fll 
(Fll+Fsl) 
Cb 

Cb 



0.97 
0.99 
0.441 
0.033 
0.406 



± 0.01 
± 0.01 
± 0.050 
± 0.007 
± 0.032 



0.95 
0.98 
0.369 
0.035 
0.406 



± 0.02 
± 0.01 
± 0.032 
± 0.007 
± 0.032 



0.95 
0.98 
0.350 
0.038 
0.413 



± 0.02 
± 0.01 
± 0.029 
± 0.007 
± 0.032 



the systematic uncertainty on this quantity. The value 
of the like-sign dimuon charge asymmetry in semi- 
leptonic 5-hadron decays is found to be 



A'' 



(-0.787 ± 0.172 (stat) ± 0.093 (syst))%. (56) 



This measurement disagrees with the prediction of the 
standard model by 3.9 standard deviations and provides 
evidence for anomalously large CP violation in semi- 
leptonic neutral B decay. The residual charge asymme- 
try of like-sign dimuon events after taking into account 
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TABLE XXIII: Values of A^^ with their statistical and sys- 
tematic uncertainties obtained for different IP selections. 



Selection 


Sample 


Central 


Uncertainty xlO^ 






value xlO^ 


statistical 


systematic 




Ifi 


-1.042 


1.304 


2.313 


All events 


2n 


-0.808 


0.202 


0.222 




comb. 


-0.787 


0.172 


0.093 




iM 


-3.244 


4.101 


7.466 


IP < 50 iJ.m 


2/i 


-2.837 


0.776 


1.221 




comb. 


-2.779 


0.674 


0.694 




1^ 


-0.171 


0.343 


0.311 


IP > 50 urn 


2/i 


-0.593 


0.257 


0.074 




comb. 


-0.533 


0.239 


0.100 




iM 


-1.293 


3.282 


5.841 


IP < 80 fim 


2(1 


-1.481 


0.541 


0.810 




comb. 


-1.521 


0.458 


0.501 




1/1 


0.388 


0.280 


0.179 


IP > 80 /jm 


2yU, 


-0.529 


0.285 


0.048 




comb. 


-0.472 


0.226 


0.091 




iM 


-1.654 


2.774 


4.962 


IP < 120 yum 


2n 


-1.175 


0.439 


0.590 




comb. 


-1.138 


0.366 


0.323 




iM 


-0.422 


0.240 


0.121 


IP > 120 


2/i 


0.818 


0.342 


0.067 




comb. 


-0.579 


0.210 


0.094 



TABLE XXIV: Mean mixing probability (xd) obtained in sim- 
ulation, and the coefficients Cd and Cs in Eq. (2), used for 
different selections. 



Sample 


Xd{MC) 


Cd 


c. 


IP<50 


0.059 ± 0.002 


0.316 ± 0.021 


0.684 ± 0.021 


IP<80 


0.069 ± 0.002 


0.351 ± 0.022 


0.649 ± 0.022 


IP<120 


0.084 ± 0.002 


0.397 ± 0.022 


0.603 ± 0.022 


IP>50 


0.264 ± 0.004 


0.674 ± 0.020 


0.326 ± 0.020 


IP>80 


0.299 ± 0.004 


0.701 ± 0.019 


0.299 ± 0.019 


/P>120 


0.342 ± 0.004 


0.728 ± 0.018 


0.272 ± 0.018 



all background sources is found to be 

Ares = (-0.246 ± 0.052 (stat) ± 0.021 (syst))%. (57) 

It differs by 4.2 standard deviations from the standard 
model prediction. 

Separation of the sample by muon impact parameter 



TABLE XXV: Measured values of afi and aj for different 
muon IP thresholds. In each column, the measurements us- 
ing the samples with muon IP larger and smaller than the 
given threshold are combined. We also give the correlation 
Pds between a^i and a^i- 



Quantity 




muon IP threshold 






50 /xm 


80 nm 


120 /um 


ai X 10^ 


+1.51 ± 0.93 


+0.42 ± 0.68 


-0.12 ± 0.52 


all X 10^ 


-4.76 ± 1.79 


-2.57 ± 1.34 


-1.81 ± 1.06 


Pds 


-0.912 


-0.857 


-0.799 



allows for separate extraction of a^j and a*). We obtain 

afi = (-0.12 ±0.52)%, 

all = (-1.81 ± 1.06)%. (58) 

The correlation pds between these two quantities is 

Pds = -0.799. (59) 

The uncertainties on a^j and a^j do not allow for the 
definitive conclusion that a^j dominates the value of A^^. 

Our results are consistent with the hypothesis that 
the anomalous like-sign dimuon charge asymmetry arises 
from semi-leptonic 6-hadron decays. The significance of 
the difference of this measurement with the SM predic- 
tion is not sufficient to claim observation of physics be- 
yond the standard model, but it has grown compared to 
our previous measurement with a smaller data sample. 
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Appendix A: Combination of two measurements 
using a scan 

In this analysis, the value of is obtained from the 
linear combination in Eq. (13). The parameter a is se- 
lected to minimize the total uncertainty on A^^, taking 
into accoimt the correlation among different contribu- 
tions to the uncertainty on A^y. This procedure is equiv- 
alent to the standard procedure of taking a weighted av- 
erage. 

To demonstrate this, we consider a model in which we 
obtain the quantity x using two measurements a and A. 
Suppose that a and A depend linearly on x: 

a = kx + b, 

A = Kx + B, (Al) 

where k, K, b, and B are parameters determined in the 
analysis, and correspond to the measurement of A^i. Us- 
ing the measurements of a and A, we obtain two esti- 
mates of x: 

Xi = {a — b)/k, 

X2 = {A-B)/K. (A2) 
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We denote the uncertainties on xi and X2 as cti and a2, 
respectively. 

Consider tlie case where the measurements of a and A, 
as well as the uncertainties cti and a2, are statistically 
independent. In this case, the value of x can be obtained 
as a weighted average: 



X = {wiXi + W2X2)/w, 

= 'i-fo-i, « = 1,2, 
w = l/al + l/a-l- 



(A3) 



Consider another estimate of x using the difference 

A' =A- aa, (A4) 

where a is a free parameter. The value of x obtained 
from Eq. (A4) is 



(A-B)- a{a - b) 
K-ak ■ 



(A5) 



Provided that the two measurements a and A, as well 
as the uncertainties ui and C72 are statistically indepen- 
dent, the minimal uncertainty on x is obtained for 



-{Kai)/{kal). 



(A6) 



The central value and uncertainty on x obtained from 
Eq. (A5), with a = amin, arc exactly the same as the 
central value and uncertainty obtained from the weighted 
average (A3). This case is similar to the combination 
(46) of the two measurements with IP > 120 iim that 
have reduced correlations. The; coefficient a in this case 
is negative, and its value depends on the uncertainties cti 
and iJ2. 

Consider another extreme case, where k = and B 
is fully correlated with &, e.g., B = Ch, where C is a 
coefficient. In this case, the value of x obtained from 
Eq. (A5) is equal to 



A- aa-{C - a)b 
K ■ 



(A7) 



Provided that (T(a) ^ o'(j4), the minimal uncertainty of 
X is obtained for amin — C. This case corresponds to the 
measurement of Eq. (36) with the full data sample. The 
value of amin is positive for C > 0. 

These two examples demonstrate that the method of 
the a scan used in this analysis is equivalent to the 

weighted average of two measurements, taking into ac- 
count the correlation among different uncertainties. 
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